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We consider a model where the strong magnetic fields of magnetars arise from 
a high baryon density, magnetized core. In this framework magnetars are 
distinguished from pulsars by their higher masses and central density. For 
magnetars, as core densities exceed a threshold, the strong interaction induces 
a phase transition to a ground state that aligns all magnetic moments. The core 
magnetic field is initially shielded by the ambient high conductivity plasma. 
With time the shielding currents dissipate transporting the core field out, first 
to the crust and then breaking through the crust to the surface of the star. 
Recent observations provide strong support for this model which accounts for 
several properties of magnetars and also enables us to identify new magnetars. 

Introduction 



Magnetars are neutron stars with surface magnetic fields ( 10 14 ( 15 ) G) . The magnetars have spin 
down ages of 10 3 — 10 5 years. Over this period, they emit a quiescent radiative luminosity of 
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10 35 - 10 36 erg/s. Besides, some of them emit repeated flares or bursts of energy typically of 
10 42 — 10 44 erg, and at times much higher (1, 2). The periods of magnetars fall in a surprisingly 
narrow window of 5-12 s. 

At such large periods, the energy emitted in both quiescent emission and flares far exceeds 
the loss in their rotational energy . The most likely energy source for these emissions is their 
magnetic energy (3, 4), yet there is no evidence of a decrease in their surface magnetic fields 
with time (5). 

There have been many attempts to explain some of this physics of which the most canonical 
is the magnetar model of Duncan and Campbell (3, 4), which is otherwise known as the dynamo 
mechanism for magnetars. This model requires the collapse of a large mass progenitor to a 
star which starts life with a period close to a millisecond. Such a fast rotation can amplify 
the inherited pulsar valued field of 10 12 G to 10 15 G. However, as described below, several 
observations on magnetars are hard to understand from such a model. 

Many authors, Leahy et al (<5), Wickramsinghe et al (7), Gaensler et al (8) and Muno (9), 
have observed that the magnetars are descendents of large mass progenitors or protostars. Wick- 
ramsinghe et al (7), based on simulations, further suggest that magnetars may actually belong 
to the higher ( than pulsars) mass (~ l.QM sun ) population of neutron stars.This in itself would 
be consistent with the dynamo model as stars descended from large mass progenitors have the 
smallest periods. However, as deduced by Vink et al (10), from the explosion energies of the 
Super Nova remnants (SNRs) of magnetars and Heger et al (11), using stellar evolution codes, 
even such massive stars will be born with periods (~ 5ms ) somewhat larger than what is 
required for dynamo amplification. 

The spin down age of magnetars deduced from their pulsar like dipole radiation is not the 
real age for the magnetars, which is instead given by the age of their supernova remnants(SNR) 
and is much larger ( see (12, 73)and (6)). Leahy et al (6) point out that this can be explained if it 
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takes an additional time of the order of 10 4 yrs for the magnetar surface magnetic field to attain 
its high final value. They interpreted it as a delayed amplification. Such a delayed amplification 
would not conform to the dynamo model. 

If magnetars are born with thier high magnetic fields, then, after a flare, one would expect 
a decrease in magnetic field and energy of the magnetar accompanied by a fall in the dipole 
radiation mediated spin down rate. However, the opposite is seen: the spin down rate increases 
after the flare indicating an increase in the surface magnetic field ( see (12-15)). Further, inspite 
of the steady X-ray emission, the magnetic field of such magnetars appears to remain high all 
the way till the end of spin down, when the stars have the largest periods. 

These striking discrepancies with the dynamo model indicate that we need an alternative 
model for magnetars. 

In an earlier work (16) it was shown that it may be possible to explain many unusual features 
of magnetars if they have a core with a large magnetic moment density, created by the strong 
interaction. Initially the core magnetic field is shielded by the electron plasma in and around 
the core. In time, the shielding currents dissipate till finally the field emerges at the surface of 
the star. 

In this work we break new ground. In section II we give a new two stage description of 
the emergence of the shielded core field from the core to the crust and then from the crust 
to the surface.This includes a consistent energy budget to set up the shielding currents and 
their subsequent decay into neutrinos, quiescent radiation and flares. Section III is new and 
devoted to detailing the extended observational support for the model. In section IV we use this 
model to establish new criteria which enable us to identify new magnetars in a list of several 
high magnetic field pulsars, whose magnetic fields are smaller than assumed for magnetars ( 
10 14 ( 15 )G) or periods much less than typical magnetars periods. 
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The Model 



Pulsars, which have radii of ~ 10km and surface magnetic fields of ~ 10 10 — 10 12 G, are 
believed to inherit such fields from their progenitors, which are stars of radius ~ 10 6 /cm and 
magnetic fields of ~ 1 — 10 2 G due to 'conservation' of magnetic flux during stellar collapse 
(17). We call such fields inherited (some authors call them fossil fields). 

Our starting point, to make the distinction between pulsars and magnetars, is that most 
observed pulsars have masses of the order of 1.4 solar masses. Theory allows pulsar(neutron 
star) masses in the larger range of 0.1 - 2.0 solar masses. Pulsars with masses larger than 1.4 
solar masses are also known (18). Magnetars are expected to be larger mass stars, whose core 
density exceeds about three times the nuclear density. This leads to a new strong interaction 
ground state -a neutral pion condensate - which aligns the spins of the neutrons (quarks) in the 
core producing a very large dynamical (spin) magnetization of the core, which can give core 
magnetic fields as large as 10 16 - 10 17 G (19, 20). 

As the strong interaction phase transition starts aligning the magnetic moments of the neu- 
trons(quarks) the electron and proton plasma in the core responds to the local change of flux, 
shielding the magnetic field in accordance with Lenz's law. This process goes on till the phase 
transition in the core is completed and all magnetic moments are aligned.The shielding cur- 
rents in and around the core confine the magnetic field inside a volume of the order of the core. 
The typical energy in the shielding currents is theN of the order of the magnetic energy of the 
core. Eventually, the Lenz currents dissipate establishing the full unshielded dipolar field in and 
outside the core . In the process of this dissipation there is a characterstic time during which 
ambipolar diffusion (21, 22) carries the core field to the crust. The magnetic field then breaks 
out of the crust to power the radiative emissions from magnetars. 

Creation Of The Core: In our model magnetars belong to the higher mass population of 
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neutron stars. They are distinguished by high density cores which go through a phase transition 
to a spin aligned ground state that carries a large magnetic moment. One likely possibility is a 
spin aligned neutral pion condensate. Such a high density ground state was first proposed by 
Dautry and Nyman (23). A simple way to understand how such a magnetic state can occur, 
for example, in a star made of neutrons is as follows. The pion condensate standing wave 
introduces an extra term — g A q ■ <?r 3 /2 in the in the strong interaction Hamiltonian ( chiral sigma 
model ) for nucleons, where q is the condensate wave vector, a the spin operator of the nucleon 
(quasiparticle) and r 3 the third component of the isospin. For neutrons the third component of 
the isospin is negative, and thus, if the spins of the neutron quasiparticles align antiparallel to 
the direction of the condensate wave vector, the energy is considerably reduced. 

Baym (24), has a simple exposition on the neutral n condensed state in the chiral sigma 
model based on the work of Dautry and Nyman (23). The ground state energy for the aligned 
state is given by summing over all (one spin) states till the fermi momentum. The minimization 
of the ground state energy with respect to the condensate wave vector, q, yields a, q, that is pro- 
portional to the baryon density. The ground state energy gets a term from the condensate that is 
negative and proportional to, n 2 , where n is the neutron number density. On the other hand, the 
kinteic energy term for the nucleons goes as, n 5 / 3 . As the density goes up the negative conden- 
sate term dominates leading to a phase transition to the neutral n pion condensed groundstate 
for large baryon density. However, inclusion of other nucleon-nucleon interactions like the ten- 
sor interaction and the hardcore may oppose such a spin alignment of the ir condensed ground 
state, whereas isobars may enhance it. Further work is needed to find a conclusive answer. 

The same phase transition to the neutral ir pion condensed ground state occurs for the 
equation of state (EOS) for two flavour quark matter (25). However, as the quarks are taken 
to be point particles, in this case we do not have tensor and hard core forces that can undo the 
spin alignment of the ir condensed ground state. This ground state has been has been analyzed 
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in (25), who use a ground state theorem to find substantial magnetization for the n condensed 
ground state ( see eqn. 50 - 52 in ref (25)). 

In previous papers (16, 19, 20) we investigated the ground state OF composite quark and 
nuclear matter stars. The equation of state (EOS) and magnetic dipole moment density for a 
neutron star with quark matter core having a neutral n pion condensed ground state is calcu- 
lated in (19). The energy per baryon for the neutral n pion condensed groundstate and the usual 
fermi sea (the uniform, q = state) is plotted in fig 1(b) in (19)). This is reproduced in Figure. 
1, from which we can read out the energy released in the core phase transition to the neutral ir 
pion condensed state to be of the order of 30 Mev (4.5 • 10~ 5 ergs) per baryon. 
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Figure 1 : The Energetics of The Neutral Pion Condensate Transition(From Bhattacharya and 
Soni (19)) 

We note that a quark matter core occurs (19) only if the nuclear matter quark matter occurs 
at low pressure otherwise the higher pressure nuclear exterior will squeeze out the softer , rela- 
tivistic quark interior. Such conditions obtain if the minimum in the energy per baryon, E BN vs 
n B , the baryon density, in the nuclear phase is close to the minimum in the energy per baryon, 
in the quark matter phase, E B q vs n B ; for then, the slope of the tangent to the two minima in 
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the, E B vs ( 1/ub) diagram ( the so called Maxwell construction ), which gives the pressure, is 
small. 

For illustration we consider a star composed entirely of neutrons with an average core den- 
sity of 5 times nuclear density (~ 10 15 gcm~ 3 ).In a unit volume there are some ~ 10 39 neutrons 
each with a magnetic moment of fi N ~ 2 • eh/2m N c ~ 1CT 23 (CGS units) giving a mag- 
netic moment density of m ~ 10 16 (CGS units). This would result in a uniform core field 
of-B = 8ir/3 ■ m ~ 5 - 10 16 G. With the reasonable value of 3 kms for the core radius, this 
would give a surface magnetic field of ~ 10 15 G. A similar result follows for a composite star 
which has a quark-matter core (16). Unlike magnetic fields whose origins are currents, this field 
built out of strong interaction mediated alignment of magnetic moments is indestructible for all 
practical purposes. 

Magnetic Field Evolution - The First Stage: The ground state of the core, above, is the 
most likely (variational) possibility without any clinching proof. For the rest of the paper we 
shall assume that the strong interaction creates a magnetized core and then discuss the evolution 
of the magnetic field from such a core by the dissipation of shielding currents. Time scales for 
these processes, relevant to a neutron, proton, electron plasma in the interior of a neutron star, 
have been worked out by Goldreich and Reisenegger (21). Their estimates show that for typical 
temperatures in the neutron star interior, the ohmic dissipation time scale would be very long: 
> 10 11 y, and hence irrelevant. Ambipolar diffusion will, however, play a very important role, 
with a dissipation time scale of 10 4 ■ ■ T a5 6 years, where B ie is the local magnetic field 
strength in units of 10 16 G and T 8 5 is the temperature in units of 10 8 5 K, a typical value in the 
interior of a very young neutron star. 

Ambipolar diffusion lets the magnetic field move out from the core to the crust. The high 
conductivity interior plasma in the star is opaque to photons and only the neutrinos are likey to 
escape (27). The temperature must be higher than ~ 10 8,5 K to have neutrino emission as the 
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dominant energy release mechanism. The photon emission which is trapped will keep heating 
the star interior to this temperature but neither the luminosity of the star or the surface magnetic 
field will change till the magnetic field crosses out of the crust to the surface. At temperatures 
> 10 8 5 K and fields of 10 16(l5 ) G in the interior the ambipolar diffusion formula gives a typ- 
ical travel time of 10 4 ^ years from core to crust. This phenomenon (absent for pulsars) will 
manifest as a delayed amplification of the magnetar activity and will provide an observable win- 
dow to the core. It is important to emphasize that ambipolar diffusion which carries the fields 
outwards is also accompanied by dissipation of magnetic energy. A fine balance is necessary 
between the rate of this dissipation and the rate of neutrino cooling to keep these temperatures 
for as long as 10 4 yrs. This would also result in a much higher surface temperature for magnetars 
compared to pulsars. 

The Second Stage: As the strong field moves through the outer crust, mechanical distur- 
bances of the crust are likely to be triggered by the magnetic pressure, leading to glitches and 
flares. The upper crust would be unable to support stresses for magnetic field difference across 
the crust of > 10 13 G, as the maximum stress that the crust can support is estimated to be ~ 10 27 
dyne cm~ 2 (28). Only after the core magnetic field penetrates the crust does the radiative emis- 
sion and serious spin down begin. The magnetars exhibit spin down ages all the way upto lO 4 ^ 
years ( a time scale similar to that for ambipolar diffusion). Over this period, the dissipation of 
the shielding currents, is responsible for their radiative luminosity and flares, till the field attains 
its final dipolar configuration. 

As per the dynamo mechanism the magnetic field of a magnetar is expected to decrease with 
time from the dissipation of magnetic energy, though, there seems to be some evidence to the 
contrary (5). Our analysis of the timing data of Livingstone,Kaspi and Gavrill (29) also supports 
this. Such a trend is also seen in certain glitching radio pulsars (30). 

Once all the screening currents dissipate away and the magnetic field reaches its fully re- 
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laxed configuration, there is no more energy available to power the magnetar activity. This 
could explain the upper cutoff in the spin period magnetars of ~ 12 s. 

The Energy Budget: From the observed energy release by magnetars we have to account 
for the steady X-ray luminosity of 10 35 ergs/sec during the spin down age (upto 10 4 ^ 5 ^ y.) and 
several flares of io 41 ~ 45 ergs. Before this phase considerable energy release is likely from 
neutrino emission (31). A total energy in excess of 10 47 ergs is released. 

Since the shielding currents are resposible for confining the dipolar magnetic field in and 
around the core, a lower bound on the amount of energy locked into the shielding currents can 
be estimated to be of the order of field energy of the core. For a core field of ~ 10 16 ( 17 ) G and 
a core of 3 km this works out 10 48 ^ 49 ^ ergs. Even if most of the energy (31) is dissipated in 
neutrino emission during the ambipolar transport phase, we are left with enough to account for 
the above processes. The energy budget works out well. 

That still leaves open the question of the source of energy required for setting up the shield- 
ing currents. The energy released in the phase transition the neutral pion condensate phase is 
~ 30 Mev (4.5 • 1CT 5 ergs) per baryon (see fig 1(b) in (19)). This energy source works out to 
an energy of ~ 10 49 ( 50 ) ergs, for a core of 3 Km of baryon density l/m _3 .This is the driver for 
making the shielding currents and this matching of energies is rather remarkable. 

Observational Support 

We now examine the broad observational support that recent data has provided for our model. 
Progenitor Mass: Since there is a one to one correlation between the mass of the star and the 
central density, in our model, only those neutron stars that exceed a certain threshold mass, M T , 
will make high density magnetic cores that can lead to magnetars. The associated progenitor 
mass also has to be high enough to yield a higher mass neutron star. The density of the magnetic 
core and hence its magnetic moment density and core magnetic field will all go up as the mass, 
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M T , goes up - till, at the maximum mass, M max the star becomes unstable (16). Many au- 
thors have observed that the magnetars are descendents of large mass progenitors or protostars. 
Evidence for this is given in (6-9). Wickramasinghe et al (7) further suggest, based on their 
simulations, that magnetars may actually belong to the higher ( than pulsars) mass (~ 1.6M ) 
population of neutron stars. We take these as indicators that the magnetars are more massive 
than pulsars. 

Delayed Amplification : Dar and De Rujula (12-14), and Marsden et al (75), find that the 
dynamo model of SGRs cannot explain their ages (15). They also find that the spin down ages 
are much smaller than the ages of their SNR's. 

The transport of the field to the surface of the star happens in two stages according to our 
model: i) the transport to the crust whose time scale is determined by ambipolar diffusion ( 
10 4 ( 5 ) years) and ii) the emergence of the field out of the crust which sets the spin down scenario 
for the magnetar. The first stage is a hidden stage and is not captured by the spin down age. 
Thus in our model there is naturally a delayed amplification of the surface field. 

The evidence given below matches exactly with the scenario sketched out above. Leahy and 
Ouyed (6) find that the spin-down ages of magnetars are systematically smaller than the ages 
of the associated super-nova remnants(SNR). While the spin-down ages of their sample varies 
between 1000-5000 yrs, the mismatch with the corresponding SNR ages is systematically of 
the order of 10 4 yrs.. This is illustrated in Fig. 1 of their paper (6) which is reproduced as Fig. 
2 below. 

The left hand panel plots the SNR ages along with a best fit curve. The right hand panel 
plots the spin-down ages also along with a best fit curve. To facilitate comparison, the best fit 
curve for the SNR ages is included in the right panel also. One sees a clear and systematic shift 
between the two best fit curves, with the SNR's systematically older. This evidence fits well 
with our picture of the delayed amplification coming from ambipolar diffusion including the 
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Figure 2: The observational evidence for delayed amplification (From Leahy et al (6)) 

time scale alluded to before. 

Anomalous Spin-down: In the conventional dynamo mechanism the star is born with a 
large surface magnetic field. Energy release at the surface then must come at the expense of 
the magnetic energy, which would produce the opposite effect: reduce the magnetic field and 
consequently the spin down rate. 

In our model, energy is released as the crust is cleaved by the emerging magnetic field, as 
the shielding currents in the crust dissipate. Subsequently the core field exits to the surface 
increasing the surface magnetic field and the spin down rate. 

This is particularly evident in a flare phenomenon, as pointed out below by Dar and DeRu- 
jula where the release of energy in the flare is accompanied by an accelerated spin down rate 
and a consequent rise in the surface magnetic field, contrary to what the dynamo model would 
predict. According to Dar and DeRujula (12, 13), in the magnetar model of SGRs the radiation- 
energy source is the magnetic field energy. The spin-down rate of SGR 1900+14 roughly dou- 
bled from around the time of its large flare on 27 August 1998 (15, 32). This sudden doubling 
of magnetic energy is what they termed energy crisis. This has also been stressed by (15). For 
AXPs too the dynamo model faces similar difficulties. 

These disparate features of magnetars and the fact that the progenitors have higher masses 
and consequently higher fields are intrinsic to our scenario. So is the delayed amplification 
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of their magnetic fields. Further, once the core field has fully emerged after dissipation of the 
shielding currents there is no further energy to release. This is the reason that magnetars have 
spin down ages between 10 4 to 10 5 years followed by cessation of activity. 

Predictions 

Extra Neutrino emission: We have pointed out that during the ambipolar diffusion stage ( 
10 4 ( 5 ) years) most of the energy dissipation takes place via neutrino emission. For magnetars 
in contrast to pulsars, this is an additional source of neutrino emission. This would be a new 
observational signature for our model. 

New magnetars from a crossover criterion: For older or larger period (more than 5s.) 
neutron stars the new feature for magnetars is that E x , the rate of steady X-ray emission must 
dominate over E R , the rate for rotational energy loss from dipole radiation. 

According to our model, the strong magnetic field emerges out of the star much after its 
birth, unlike in the case of ordinary pulsars. We shall give some estimates which, though quali- 
tative, capture some of the physics. In our picture, when the shielding currents eventually reach 
the crust, the core magnetic field is completely shielded within the crust. As the shielding cur- 
rents in the crust dissipate, they restore the crust magnetic field to what its unshielded value 
would have been. 

Since this field was confined by the Lenz currents in the crust, the lower bound on the field 
energy carried by the shielding currents in the crust is given roughly by, E stored ~= Bj/8n ■ 
AttR 2 ■ AR cr where Bf is the final observed surface field, R the radius of the star and AR cr 
the thickness of the crust. This energy will be dissipated and released as X-rays, flares etc. in a 
typical time-scale of the spin-down age t sd . Thus, an estimate for the average X-ray luminosity 
is given by 

Lf c = k ■ R 2 AR cr -Bj-^ 
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where we have introduced a factor k as an indicator of possible deviations from this estimate., 
which is a lower bound. 

As a matter of fact we can estimate the ratio of the calculated X-ray luminosity L c x lc to E R . 
From our earlier discussion one sees that both these are proportional to P and hence their ratio 

r calc h, 

E R 4000 13 

is independent of the mechanism for spin-down. In our model both B and P go up with time. 
Therefore 1Z will keep increasing with time giving rise to a crossover between 1Z < 1 and 
1Z > 1 regimes. If we use k ~ 20 as preferred by observed data, then this ratio exceeds unity 
even for B± 3 ~ 3 and P ~ 5s. For pulsars L x is absent. 

In table I we have summarized the observed properties of a number of high magnetic field 
pulsars. The magnetic fields in table are inferred from the observed P, P by taking k — 1 and by 
the use of the dipole radiation formula . According to the crossover criterion those stars marked 
as M in table I are candidates for magnetars in our scenario. This crossover happens when the 
rotational energy loss rate (at P ~ 5s.) is too small to fuel the X-ray flux. We therefore need 
another energy source for the X-radiation - in our model it is the dissipation of the shielding 
currents. Pulsars do not have this source. This crossover during the evolution of a star is a 
signal for magnetars. We emphasize that in the conventional dynamo scenario they would not 
be considered as magnetars. We note that the object J 1846-0258 is anomalous in many respects 
and needs special consideration ( see below). 

New Young (sub second period) magnetars: 

All 'observed' magnetars have periods between 5 < P < 11 sec, magnetic fields between 
10 14 < B < 10 15 G and a steady radiative emission of ~ 10 35 ergs/sec. Why, then, are 
magnetars not observed with P < 5s. ? 

In the dynamo model magnetars are born with high B fields ~ 10 14 ( 15 ) G and would spin 
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down very fast from periods of a few ms. to periods of over a second(in a few hundred years). 
Therefore there is a much higher probability of observing them when they have slowed to pe- 
riods ~> 5 s. and their spin down ages are much higher. In this epoch their X-ray luminosity 
L°x s is much greater than Er. Hence it is customary to think of magnetars in the above regime. 

Observation of PSR 1846-0258 with P=0.326 s., B = 5 • 10 13 G and L$ s < E R would 
then rule against it being identified as a magnetar conventionally. Furthermore, since the dipole 
radiation dominates over the staedy X-ray flux, the crossover criterion discussed earlier can not 
be applied to it at this epoch. 

The pulsar PSR 1846-0258 has a very high X-ray luminosity, comparable in magnitude 
to that for magnetars. One may then suspect that it may become a magnetar when its period 
evolves to P ~ 5s. We need a new signature to establish its credentials as a magnetar. 

• Gavrill et al (33). report a glitch accompanied by a flare in this object that are more 
characterstic of magnetars than pulsars. The post flare phenomenon shows an increase in 
the spin down rate. 

• A careful and detailed timing data and analysis for this pulsar has been done by Kaspi et 
al (29). Part of their analysis most relevant to our model is best summarized in Fig. 4 
of (29) which we have reproduced in Fig. 3 here. 
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Figure 3: Timing Analysis from Kaspi et al (29) 
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Prior to the glitch the frequency derivative z> was —6.3 ■ 10~ n s -2 . Immediately after the 
glitch this changed to —7.55 ■ 10 _11 s -2 . They observed an exponential decay of this increase 
with a characterstic time scale r d = V27d. But even after 920 d.(7.25r d ) they observed an 
asymptotic v of —6.65 • lCT n s~ 2 giving a difference of 3.5 ■ lCT 12 s -2 over the preglitch value. 
The exponential tail of the post-glitch increase after 7.25^ being 5.5 • 10~ 14 s~ 2 , we interpret 
their data as evidence for a real change in v of 3.5- 10~ 12 s~ 2 . Using the relation B = 3.2y(P P- 
10 19 G which is based on the dipole formula 

• _ 2 4tt 2 B 2 R 6 
~ 3~P c 3 I 

we can translate this asymptotic change in v to an increase in the surface magnetic field of 
0.13 ■ 10 13 Gin about 2.5 years. We have made use of the fact that P only changed by a few 
ppm during this entire period. 

A post flare anomalous increment in the surface B and the spin-down rate, is predicted 
by our model for magnetars. In the context of our model the anomalous behaviour of PSR 
1846-0258 identifies it as a magnetar. 

It is important to see how distinctive PSR 1846-0258 is in this regard when compared to 
other glitching pulsars which are obviously not magnetars. In particular, whether they too 
have persistent changes in v that are significant, in which case persistent changes in v would 
merely reflect the large glitches that took place and nothing more. We have analysed the cases 
of the Vela pulsar and PSR 0355+54 both of which are known to have large glitches with 
Af2/f2 ~ 1CT 6 . A comparison of these two cases with that of J1846-0258 is given in Table. 3. 
In the case of PSR 0355+54 (34), Lyne's fit (35) for the observed post-glitch behaviour yields a 
persistent shift Atl/tl = 0.0059 which is ten times smaller than in the case of PSR 1846-0258. 
In this case observations were carried out over a 77 d. period while the characterstic timescale 
of the exponential decay was 44 d. In the timing analysis of Kaspi et al (29) the shift was 
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directly observable as data was available over a much longer period. In the case of the Vela 
pulsar, according to Alpar (36), the observed long term breaking index of 1 .4 can be explained 
by a persistent shift in Atl/Cl of 3 x 1CT 4 happens about once every four years. Even this 
inferred effect is more than two orders of magnitude smaller than the observed effect in PSR 
1846-0258. This indeed makes PSR 1846-0258 distinctive. 

New braking index for J1846-0258: In a very interesting development Kaspi et al (57) 
report a new braking index of 2.16± .13 for the post-outburst phase of this object which is 
considerably lower than the pre-outburst value of 2.65±.01. They carefully analyse the sources 
for such a lowering of the braking index and essentially identify three such: 

• A positive second time-derivative of the moment of inertia I. They conclude that it is hard 
to imagine a physical situation causing this. 

• A positive rate of change of a, the angle between the magnetic axis and the axis of rota- 
tion. However, according to them no change in the pulse profile has been observed over 
the relevant period, making this source very unlikely. 

• Plasma distortions of magnetic fields can also be a source. But according to these authors 
that would result in a six fold increase in particle luminosity and no such increase has 
been observed. 

• An increase in the surface magnetic field B, more precisely a positive 

Thus a persistent increase in the surface magnetic field is the most likely explanation for the 
lowering of the braking index and this further corroborates our point of view that the surface 
magnetic field of J1846-0258 is indeed increasing and that this object is a magnetar in making. 

Another possibile source for an increase in v is pulsar wind - but there is no reason to expect 
it to correlate with long term post flare phenomena. 
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Conclusion 



In conclusion we enumerate some of the consequences of the model presented above: 

i) Magnetars belong exclusively to the higher than pulsar mass population of neutron stars 
that have a high density magnetic core. 

ii) This core is created by the strong interactions and the most likely ground state is a ir 
condensate that aligns magnetic moments to create large dynamical 10 16 ( 17 ) G magnetic fields 
at the surface of the core. Dynamical fields are 'permanent' unlike fields derived from currents. 

iii) The core field is shielded by Lenz currents generated in the high conductivity plasma in 
and around it, but is gradually transported to the crust by ambipolar diffusion over a timescale 
of ~ 10 4 years - this results in a time delay before the field comes out to the surface. 

iv) The strong magnetic field breaks through the crust as the shielding currents dissipate 
giving out a steady X-ray flux and several energetic flares. 

v) This further implies that the surface field keeps increasing in magnitude till all shielding 
currents dissipate and the permanent dipolar core field is established. 

We have found that all these phenomena are supported by extensive data and observations. 
We have also used the model to identify some stars, in a sampling of a few high magmetic field 
pulsars, as magnetars. Thus this model for magnetars throws up a lot of unexplored physics 
from the strongly interacting core to the plasma physics and the crustal solid state physics of 
huge magnetic fields. 

Our understanding of neutron stars is at a crossroad. We have to understand many families 
of neutron stars, for example pulsars and magnetars, in one framework. This is what we have 
tried to do in this work. Neutron stars are also the laboratory to understand the high density 
phase diagram for strong interactions. This work gives us a new understanding of the strong 
interactions that is linked intimately to astrophysical data. 



17 



Acknowledgements 



V.S. thanks Dipankar Bhattacharya for earlier collaboration that led to this work . We are 
grateful to A. Deshpande and Raymond F. Sawyer for a critical reading of the manuscript and 
for insightful comments that have significantly improved the paper. We thank S. Jhingan, Pankaj 
Sharan, C. Sivaram, S. Mahajan, B. Paul, Deepak Kumar, P. Dasgupta, A. Khare, G. Baskaran 
and D. Sahdev for discussions. Special thanks are to Siraj Hasan for facilitating initial work at 
the Indian Institute of Astrophysics. V.S. thanks the University Grants Commission, National 
Physical Laboratory and Centre for Theoretical Physics, Jamia Milia University for support and 
the CHEP, IISc, Raman Institute and Indian Institute of Astrophysics for hospitality during the 
course of this work. NDH wishes to acknowledge the Department of Atomic Energy, India, for 
the Raja Ramanna Fellowship and CHEP, IISc, for the opportunity to utilise this fellowship. 

Note Added 

After this work was completed and this paper written we came across a paper by Kutschera (39) 
which has also discussed shielding and eventual outward transport of core magnetic fields. But 
the mechanism discussed there is very different from Ambipolar diffusion, and that paper is 
addressing low field millisecond pulsars. It estimates much larger time scales and will not be 
relevant for magnetars. 
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TABLES 



Table 1: Compilation of actual data, P, P, for high magnetic field pulsars (HMFP) from ref. 
(38), with inferred B, tsd, Er ( from the dipole radiation formula). L^' c ( = E stor ed/TsD) is 
calculated for k = 1. In estimating Lx we have taken a star radius of 10 kms, a crust thickness 
of 1 km and a moment of inertia I of ~ 10 45 gem 2 . Our estimates for L x are about an order of 
magnitude smaller than the observed values. This means that the A;-factor is in reality about 20. 



Name 


P 

in sec. 


Pin 

io- 12 


t S d in 
10 3 yrs 


B 13 


E R in 
10 33 ergs/s 


Lf s in 
10 33 ergs/s 


Lf lc in 
10 33 ergs/s 


J^obs I j^calc 
X 1 X 


J1718-3718 M 


3.3 


1.5 


34 


7.4 


1.6 


6 


0.25 


24 


Jl 846-0258 M 


0.33 


7 


.7 


5 


8300 


100 


5 


18 


Jl 119-6127 


0.407 


4.1 


1.6 


4.1 


2600 


0.216 


1.7 


0.13 


Jl 847-0130 M 


6.7 


1.3 


83 


9.4 


0.17 


3 to 8 


0.2 


19 


J1814-1744M 


4.0 


0.74 


85 


5.5 


0.47 


4.3 


0.06 


72 


PSR0154-61 


2.35 


0.189 


197 


2.1 


0.57 


0.14 


0.004 


35 
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Table 2: A comparison between pulsars at the beginning of spin-down, high magnetic field 
pulsars (HMFP)and typical magnetars. Other parameters are calculated using the dipole formula 
with P, B as input.In the first row we have given the typical values of P, B at the beginning of 
the spin down phase.The second row is constructed by taking the average values of P, B from 
table I (HMFP)for all cases marked M (except J1846-0258). Finally, the last row is obtained by 
taking for P, B from the average of the AXP data given in Gonzalez et al (38). 



Name 


P 

in sec. 


P 

in 10- 12 


T SD = P/2P 
in 10 3 yrs. 


Er 

in ergs. 


Er 

in ergs/s 


-Bl3 


j^calc 

in ergs/s 


Beginning of spin down 


0.1 


1.0 


1.666 


2 • 10 48 


4 • 10 37 


1.0 


10 32 


Typical High Mag Field HMF 
pulsar or magnetar: P > 5s 


5.0 


1.13 


74.0 


8 • 10 44 


3.6 • 10 32 


7.5 


10 32 


Typical magnetar P > 5s 


9.0 


2.5 


60 


2.5 • 10 44 


1.25 • 10 33 


50.0 


6 • 10 34 



Table 3: Comparison of persistent A(l/(l of PSR 1846-0528 with those in some large-glitch 



Name 


Length of 
observation 


Exponential 
decay time Td 


Persistent 

An/n 


Source 


PSR 1846-0528 


920 d. 


127 d. 


5.5- KT 2 


Timing Data (29) 


PSR 0355+54 


77 d. 


44 d. 


6 • 10" 3 


Timing Data (35) 


Vela 






3 • 10" 4 


Breaking Index (36) 
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